Test protocole d’étude des contrainte et objectifs après scénarisation
Clémentine de Montgolfier
2025-08-19
## # A tibble: 277 × 4
## folder plot_id simul simulation_id
## <chr> <chr> <chr> <int>
## 1 output-cmd_1.txt RETZ_00972_01 simulation_10 1
## 2 output-cmd_1.txt RETZ_00972_01 simulation_11 2
## 3 output-cmd_1.txt RETZ_00972_01 simulation_12 3
## 4 output-cmd_1.txt RETZ_00972_01 simulation_13 4
## 5 output-cmd_1.txt RETZ_00972_01 simulation_14 5
## 6 output-cmd_1.txt RETZ_00972_01 simulation_15 6
## 7 output-cmd_1.txt RETZ_00972_01 simulation_16 7
## 8 output-cmd_1.txt RETZ_00972_01 simulation_17 8
## 9 output-cmd_1.txt RETZ_00972_01 simulation_18 9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19 10
## # ℹ 267 more rows## # A tibble: 278 × 4
## folder plot_id simul simulation_id
## <chr> <chr> <chr> <int>
## 1 output-cmd_1.txt RETZ_00972_01 simulation_10 1
## 2 output-cmd_1.txt RETZ_00972_01 simulation_11 2
## 3 output-cmd_1.txt RETZ_00972_01 simulation_12 3
## 4 output-cmd_1.txt RETZ_00972_01 simulation_13 4
## 5 output-cmd_1.txt RETZ_00972_01 simulation_14 5
## 6 output-cmd_1.txt RETZ_00972_01 simulation_15 6
## 7 output-cmd_1.txt RETZ_00972_01 simulation_16 7
## 8 output-cmd_1.txt RETZ_00972_01 simulation_17 8
## 9 output-cmd_1.txt RETZ_00972_01 simulation_18 9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19 10
## # ℹ 268 more rowscomplete_scene <- import_output_scene( output_name = “complete”, base_path = “C:/Capsis4/data/forceps/clementine/Test_protocole/”, output_file = “output/protocole_completeScene.RData”) %>% filter(date <= “2100”) %>% group_by(plot_id, date, scenario) %>% summarise( totalVolumem3.ha = sum(totalVolumem3.ha, na.rm = TRUE), totalBiomass_t_ha = sum(totalBiomass_t_ha, na.rm = TRUE), .groups = “drop” )
Etats initiaux des placettes
| Plot ID | Surface | Type de peuplement | Essence déterminante | Essence dominante 1 | Essence dominante 2 | Structure et occupation du sol | Surface terrière (m2/ha) | Diamètre moyen (cm) | Âge médian (ans) |
|---|---|---|---|---|---|---|---|---|---|
| RETZ_00230_03 | 2.50 | FHETE | HET | HET | F | 0.0 | 7.5 | 12.0 | |
| RETZ_00816_02 | 0.94 | FHET1 | HET | HET | F | 17.5 | 12.5 | 12.0 | |
| RETZ_00901_01 | 10.23 | FHETG | HET | HET | CHP | F | 21.0 | 57.5 | 144.5 |
| RETZ_00972_01 | 4.59 | FHETG | HET | HET | CHS | F | 31.0 | 57.5 | 144.5 |
| RETZ_01008_05 | 3.38 | FHETS | HET | HET | CHS | F | 0.0 | 3.5 | 12.0 |
| RETZ_01067_02 | 1.73 | FCHPE | CHP | CHP | F | 0.0 | 7.5 | 12.0 | |
| RETZ_01210_04 | 3.10 | FCHSE | CHS | CHS | F | 0.0 | 7.5 | 12.0 | |
| RETZ_01307_02 | 1.28 | FDOUG- | DOU | DOU | F | 29.0 | 52.5 | 59.5 | |
| RETZ_01368_01 | 7.91 | FHETM | HET | HET | CHA | F | 31.0 | 37.5 | 114.5 |
| RETZ_01415_02 | 2.23 | FCHPG+ | CHP | CHP | CHA | F | 17.5 | 62.5 | NA |
Exploration des itinéraires (pour validation des dynamiques forestières)
Pour rappel des itinéraires :
- Coupe rase : (simulation1)
- coupe d’aménagement tous les 10 ans 30% (type de coupe 0.5) en gardant la composition initiale
- coupe rase à l’âge de maturité de l’essence majoritaire
- Irrégulier mixte : (simulation2)
- coupe partielle tous les 5 ans de 10% (type de coupe 0.5)
- en cherchant une composition égale des 4 espèces principales
- Sans coupe : (simulation3)
- libre évolution sans intervention
Je ne suis pas sûre de pourquoi il y a des chutes de surface terrière dans la simulaion 3 (sans coupe), l’itinéraire étant bien : 150_3_0.5_0%_FSyl-80. peut-être juste mortalité naturelle
Calcul et visualisation des indicateurs par placette
Sur les 80 ans de simulations on calcul les indicateurs suivants :
Indices de diversité de Hill (richesse, Shannon, Simpson, Hill 0, Hill 1, Hill 2)
Volume extrait par scénario
Biomasse sur pied totale et par espèce
Indicateurs globaux
Relations entre indicateurs globaux
